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Abstract. Evolutionary synthesis models (ESM) have been extensively used to obtain the star 
formation history in galaxies by means of SED fitting. Implicit in this use of ESM is that (a) for 
given evolutionary parameters, the shape of the SED is fixed whatever the size of the observed 
cluster (b) all regions of the observed SED have the same weight in the fit. However, Nature 
does not follow these two assumptions, as is implied by the existence of Surface Brightness 
Fluctuations in galaxies and as can be shown by simple logical arguments. 
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^ . 1. The Deterministic ESM problem 

I \ Evolutionary synthesis models (ESM) have been extensively used to obtain the star 
^ ' formation history in galaxies by means of Spectral Energy Distribution (SED) fitting 
CNJ techniques or the position of galaxies in theoretical diagnostic diagrams. This statement 
is confirmed by most of the contributions in these proceedings, and by the scientific liter- 
ature. Implicit in this use of ESM here is a deterministic interpretation of model results: 
the shape of the SED is univocally determined by evolutionary conditions (birthrate, 
J — metallicity and age), and does not depend on the size of the system and the spatial 
resolution of the observation. 

There are, however, strong evidences that Nature does not obey this simple assumption. 
Evidences range from simple logical reasoning to observational facts, as the existence of 
' Surface Brightness Fluctuations in distant galaxies. 

^ The simple logical reasoning can be formulated in the following terms. For the sake of 

' reasoning, let us assume a single stellar population (SSP) model (i.e. the star formation 
Ci • rate is a Dirac's Delta function): 

^ I (a) Synthesis models can be scaled to clusters of any size, so, in principle, this model 
• ■ should also be applicable to any individual star. However, determining the cluster's evolu- 
K% ' tionary parameters by SED fitting will only produce a realistic solution in the wavelength 
range (if any) where the SSP model is dominated by stars with the same spectral type 

" " as the chosen individual star. 

(b) In the case of clusters, similarly, the realism of the parameters inferred from SED 
fitting in each particular stellar cluster depends on how Nature has sampled the different 
stellar evolutionary phases. 

(c) Finally, not all wavelengths/indices can be fitted with the same accuracy. SED 
ranges, line profiles and indices dominated by very specific stellar types only provide 
information on how many stars of the given type are in the cluster, but not on the 
overall composition of the cluster. As an example, infrared emission mainly provides 
information on stars in the giant and asymptotic giant branch: very luminous but very 
sparse (certainly sparser than Main Sequence stars) stellar types. 
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Figure 1. Left: Evolution of the moments of tlie distribution that describe the possible values 
of different photometric bands. Right: Detail of the variation of the moments as a function of 
the wavelength for a IGa old cluster. 

2. The Probabilistic ESM solution 

When ESM are placed in a statistical framework, the previous problems do not apply: 
standard ESM results only provide the mean value of the distribution that describes 
the universe of all possible values of the observable (SED, spectral index, etc.) taking 
naturally into account the size (number of stars) of the system (the population Luminos- 
ity Distribution Function, pLDF). Additionally, high-order moments of this probability 
distribution and their scale relations with the size of the system can be easily obtained. 
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As en example, we show the evolution of the mean, variance, skewness and kurtosis 
of different photometric bands as a function of age in the left panel of figure [TJ It can 
be seen that the IR bands, having a larger variance, skewness and kurtosis (i.e. more 
asymmetric pLDF) than the blue bands, are less robust than the blue bands for fitting 
the results of synthesis models with observed data. 

In the right panel we show a portion of the SED for a IGa old cluster. It can be seen 
that the line profile (A ~ 4102A), having a larger skewness and kurtosis, is a less 
robust index than the continuum or than the H and K Ca II (plus Us) profiles. 

The probabilistic framework provides a natural test to evaluate the quality of spectral 
fitting, and to avoid overfitted results. In this framework, the fitting of observational data 
would produce less precise, but much more accurate (physical) results. 

Acknowledgements 

This work was supported by the Spanish project PNAYA2004-02703. MC is supported 
by a Ramon y Cajal fellowship. VL is supported by a CSIC-I3P fellowship. 



References 

Cervifio, M. & Luridiana, V. 2006, A&A 451, 475 



